Decay of massive scalar field in a Schwarzschild background 
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£NJ ■ Abstract 

■ The decay of massive scalar field in the Schwarzschild black hole background is investigated here 
by consideration its quasinormal spectrum. It has been proved that the so-called quasi — resonant 
modes, which are arbitrary long living (purely real) modes, can exist only if the effective potential is 

■ not zero at least at one of the boundaries of the _R-region. We have observed that the quasinormal 
spectrum exists for all field masses and proved both analytically and numerically that when n — > oo 

£SJ , the real part of the frequencies approaches the same asymptotical value (In 3/(8ttM)) as in the case 

^ ■ of the massless field. 

in 

At late times, the decay of fields in a black hole background is dominated by some resonant 
quasinormal modes (QNMs) pQ. The interest in QNMs originated from possibility of detection of grav- 
' itational waves from colliding black holes, especially after creating the new generation of gravitational 

antennas. Recently the study of the quasinormal modes has gained a strong motivation coming from 
""q ' ADS/CFT correspondence, where the QNMs can be interpreted as the poles of the temperature Green 

function in the dual conformal field theory. It is also possible, that QNMs play an important role in 
Loop Quantum Gravity. All this motivated the extensive research of QNMs for different black holes and 
different fields (both massive and massless) around them (see for instance [2] and references therein). 

Long time ago, when considering the problem of a massive scalar field around Schwarzschild black 
hole, A. Starobinskii and I. Novikov found, that, when working in the in the frequency domain, the poles 
of the Green function in the complex plane are closer to the real axis than in the massless case. Thus, it 
was observed that the "massive" modes decay more slowly than the massless ones. Later it was confirmed 
in a series of papers, both in frequency domain 0, 0], and in time domain [H]-[E], that the greater the 
mass of the field, the slower the decay. Then, a natural question is, what will happen at further increasing 
of the mass of the field? Can the decreasing decay rate approach zero leading to existence of arbitrary 
long living modes? Which is high overtone behavior in this case? 

In the present paper we are answering these questions by thorough investigation of the quasinormal 
behavior corresponding to the decay of a massive scalar field. 
1. Basic equations 

The massive scalar field in a curved background is governed by the Klein-Gordon equation: 
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□$ - m 2 $ = -j== ) ,„ -m 2 $ = 0. (1) 

After the separation of angular and time variables, the result radial equation for the Schwarzschild 
background 

rlr 2 9 M 

ds 2 =f(r)dt 2 --^-r*(d9 2 - S in 2 6d<P% /(r) = 1 - ±f-, (2) 
can be reduced to the following wave-like equation 

d 2 



dr* 



+ co 2 -V(r*) Ur(r*) = 0, (3) 
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where 



and / = 0, 1, 2, 3 . . . parameterizes the field angular momentum. 

The effective potential has the form of the potential barrier which approaches some constant values 
both at the event horizon and at spatial infinity. Therefore the QNMs boundary conditions (within the 
dominant asymptotic order) have the form: 

*(r*) ~ C± exp(±iK±r*), r* — > ±oo; (4) 

Under these boundary conditions the QNMs of massive scalar field were studied in several papers in the 
WKB approximations E9- Yet, if we need to explore accurate Leaver method ^3], we have to 

deal with irregular singular point at infinity. This implies that one needs to take into consideration the 
sub-dominant asymptotic term at infinity: 



*(r*) ~ C + e lxr 'r {lMm /x \ (r,r* +oo); x = Vu 2 -m 2 . (5) 

(The sign of \ is to be chosen to stay in the same complex surface quadrant as ui.) 

Within the Leaver method we can eliminate the singular factor of the solution that satisfies in- 
going wave boundary condition at the horizon and © at infinity, and expand the remaining part into the 
Frobenius series that are convergent in the i?-region (between the event horizon and the infinity). The 
appropriate series are: 

$(r) = eixr r (2iM x+ m m s/x) ^ _ 2 j£\ £ Qn ^ _ , ( 6 ) 

\ / n \ / 

Substituting © into @ we obtain a three- term recurrent relation for the coefficients a n : 

a ai + p ao = 0; (7) 
a n a n+ i + P n a n + 7„a„_i =0, n > 0, 



where 



= (n + l)(n + 1 - 4Miu>); 



A, = M{UJ + * )(4M( " + X)2 + t{2n + + 3X)) - 2n(n + 1) -1-1(1+ 1); (8) 

X 

In = {n- Mi(uj + x) 2 /x) 2 ■ 

Since the series are convergent at infinity, the ratio of the series coefficients is finite and can be found 
in two ways: 



a n p n _i /3 „_ 2 _ a „_ 37 „_ 2 /... «n Pn+1 /3„ + 2 -a„ +27 „ +3 /. 



(9) 



Thus we have equation with respect to the eigenvalue u 

Un-lln _ Otnln+1 

a n -2ln-l ~ Q a n+l 

f) n -2— a„_37n-2/... l J n+l /3„ +2 -q„ +2 7„ +3 / 



n "-Tt— X ru "--ft rft~t--l /ia\ 

Pn ~ ~ a n -27re-i — a a„ + H„ +2 > 

P™-1 37. o-r>„ 777 777: Pn+1 



that can be solved numerically. 

To improve the convergence of the continued fraction on the right side of l|10|) we used the technique 
developed by Nollert He considered the recurrence relation 



IN _ 7jV 

Pn — ~a a ^"* N+ ~ — 7 — Pn ~ olnRn+i 



r> _ IN _ IN 

Un ~ 75 awTw+i — 75 7 — d • \ X1 ) 
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Making use of the recurrence relation Ijllfl one can find for large N: 

R N = C + dN- 1/2 + C2N- 1 + ... (12) 

where 

C = -1, 

C\ = ±2y/-iM x , Re(Ci) > 0, 
3 

C 2 = - + AMix + iMm 2 /x, 

etc. 

The series (|12[) converge for |x|/./V < A < 00, therefore we can use this approximation for Rn inside 
the continued fraction for some N 3> — Im(x) ~ n. In practice, to find an appropriate N we should 
increase it until the result of the continued fraction calculation does not change. 

2 . Quasi-resonances 

Recently Ohashi and Sakagami J2| studied QNMs for the decay of the massive scalar field and found 
that there are perturbations with arbitrary long life when the field mass has special values. They called 
these modes quasi — resonances. 

The remarkable fact is that in the massive case the purely real frequencies are not forbidden. Indeed, 
by multiplying © by ^* and integrating we obtain: 

+R2 +-R2 +R2 

J + (cu 2 - V)V) dr* = to 2 J \V\ 2 dr* - J (|*'| 2 + V\^\ 2 )dr* + = °- ( 13 ) 

-Hi -R t -R t 

Substituting the boundary conditions for large R\^, we find that 

***i:;:!J=iexp(-2Im(x) J R2)^ MIm(x)ro2/lx|2 (x\C + \ 2 +«(^)^) + (14) 

+iexp(-2hn(v)R 1 ) ( uu\C-\ 2 + o ( — ) i? d 



k -Ki- 
lt is easy to see that for the massless case (x = and lm(u) — 0, the equation l(T^|) implies 

w|C+| 2 =0 and uj\C-\ 2 = 0. 

This means that there is no real frequency in the QNM spectrum of the massless scalar field. 

When the field mass is not zero, the purely real frequencies w are not forbidden because (|13f) is satisfied 
when 

Re(x)|C+| 2 = 0, and uj\C-\ 2 =0. (15) 

There is no wave falling on the horizon (C_ = 0) in the latter case, and there is no energy transmission to 
infinity (Re(x) = 0). Therefore oscillations do not decay. The situation is similar to the standing waves 
on a fixed string. The requirement for \ to be imaginary, bounds the quasi-resonance frequencies by the 
field mass 

\uqrm\ < m. 

Eq. (|15|l also illustrates that phenomenon of quasi-resonance exists since the potential is not zero at 
spatial infinity. If the potential is zero for r* — > ±00 no quasi-resonant oscillations can exist. Thus, for 
instance, in the spectrum of the Schwarzschild-de Sitter background |17j . no quasi-resonant oscillations 
can exist even for massive scalar field. 

The numerical investigation shows, that increasing of the field mass gives rise to decreasing of the 
imaginary part of the QNM until reaching the vanishing damping rate. When some threshold values 
of m are exceeded, the particular QN modes disappear (see figures Q] and EJ) ■ We can see, that the 
larger field mass is, the more first overtones share this destiny. Thus disappearing of the QNM does not 
induce the disappearing of the whole QN spectrum. Since, as it was proved in |16j . the number of QN 
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Figure 1: First three overtones of I = QNMs of the scalar field of different field masses (M=l) 
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Figure 2: First three overtones of I = 1 QNMs of the scalar field of different field masses (M=l) 
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Figure 3: The real part of high overtones of I = 0, m = 0.3 and m = 3 (M=l) asymptotically approaches 

In 3 

modes is infinite for each multi-pole number, the presence of any finite mass of the field will not induce 
disappearing of the whole spectrum, leaving a countable number of QN modes. 

Note also, that the presence of the mass term does not change the asymptotic crucially and, in fact, 
QN frequencies found under boundary conditions Q with the 6th order WKB approximation |1(J| do 
not differ from those found under the Leaver method where sub-dominant asymptotic term is taken into 
account, at least for not large values of m (for more references on the usage of WKB method for finding 
QNMs see for instance [TT]1. 

3. High overtones 

First of all, let us propose some illustrative estimations for high-overtone behavior of the massive field. 
Starting from expansion of x a t high overtones in the following form 

and substituting (|16ll in (J8J, one can find 

a n (m) = a n (m = 0); (3 n (m) = (3 n (m = 0) - 16M 2 ™ 2 + O Q^j ; j n (m) = -y n (m = 0) + O 

We see that for high overtones (Mlm(u) ^ 1), the behavior of coefficients is the same as in massless 
case with some shifted multi-pole number I. Reminding that asymptotical behavior of QNMs does not 
depend on the multi-pole number we can conclude that 

W (m, AT) = W (m = 0, N) + O (~) = + ^ + 0(A^ 1/2 ). (17) 

If we do not pay attention to terms like n/w, this arguing illustrates why the asymptotical behavior of 
QNMs does not depend on the field mass. To prove this fact exactly, one should to reproduce all the 
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steps in |15| . We omit them because they are almost exactly the same as in |15| . but require additional 
discussions. Finally, we find the equation for high overtones that leads to the solution (|17fl : 

tan (2Miiroj) tan 2 (Miix{u + xf/x) = tan (2Miiroj) tan 2 (^IMiiruj + O (—3 J ^ =±i + 

Numerical calculations confirm our arguments. On figure [3] you can see that the difference between 
modes for m = 0.3 and m = 3 is small for high overtones. More accurately, we constructed fits for both 
curves and found that they tend to l|17|) . (Fits and data for different I are available from the second 
author upon request). 

Conclusion 

We have made a thorough investigation of the quasinormal spectrum for massive scalar field in a 
Schwarzschild black hole background. It is shown that the mass of the field has crucial influence on 
damping rate of the QNMs. In particular, the greater the mass of the field is, the less the damping rate. 
As a result, purely real modes appear which corresponds to non-damping oscillations, and, when the field 
mass is greater than some threshold values, the lower overtones disappear. This happens however with 
lower overtones only, while all the remaining higher overtones are still damping. Generally, we proved, 
that such arbitrary long living modes are forbidden, unless the effective potential is non-zero at least 
at some of the two boundaries (event horizon or spatial infinity). Note, that since we are working in a 
linear approximation and also since the QN modes do not form a complete set, the existence of arbitrary 
long living modes, apparently, does not mean we can really have a system without damping. We should 
expect that in the fully non-linear analysis the quasi-resonance modes may have non-vanishing damping 
rate. 

In addition we proved that at asymptotically high overtones the real part of the QNMs goes to 
ln3/(87rM) for any (finite) mass of the field. The numerical, analytical and WKB methods we used show 
excellent agreement in their range of validity. It is interesting for us to find out what will change in the 
above picture if considering massive fields of higher spin. The massive Dirac field represents apparently 
the most easy case. 
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